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I. INTRODUCTION 



Recently, the non-Fermi liquid infrared fixed point of the multichannel Kondo model, 
which describes a system of k identical conduction bands interacting with an impurity spin 
S, has received considerable attention [?,?,?,?,?,?,?,?,?,?,?,?,?,?,?,?,?,?]. The interest is 
aroused by its potential application in the metallic glasses [?,?,?], and heavy fermion Ura- 
nium alloys [?,?,?] where more examples have been found with their low temperature 
behavior falling out of the usual Fermi liquid expectation [?,?]. It is enhanced by the more 
interesting observation of non-Fermi liquid behavior in the normal state of the high Tc 
cuprates [?,?] and its resemblance to the low energy behavior of the multichannel Kondo 
model [?,?]. 

Although this model is more than ten years old [?] and has been attacked by various 
methods [?,?,?,?,?,?,?], there is still room left for a simple interpretation(or reinterpretation) 
of the physics of the low energy fixed point which is often said to be nontrivial. This task 
is easily achieved in the limit k ^> 1 where the low energy fixed point has a value of order 
1/k for the coupling constant and is accessible to a renormalization group(RG) improved 
perturbative expansion in 1/k. Some physical quantities difficult to calculate by the other 
methods are readily obtained by this perturbative approach and the results provide new 
insight into the fixed point. That the nature of the fixed point remains qualitatively the 
same when continuing k to small values as long as k > 2S has been demonstrated by the 
results of Bethe Ansatz or conformal field methods [?,?,?]. In this paper, we calculate a 
long list of physical quantities to the leading or sub- leading order in 1/k. A simple physical 
picture can be sketched based on these results and the previous understanding. 

For an antiferromagnetic Kondo interaction between the the impurity spin and the con- 
duction electrons, the impurity spin pulls in conduction electrons with opposite spin. This 
increases the spin-flip exchange which in turn enhances the attraction of the conduction elec- 
trons with opposite spin to the impurity spin. The result of this cooperative enhancement of 
the Kondo interaction is that one conduction electron from each of the k channels is pulled 
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in to screen the impurity spin. Since k > 2S, the conduction electrons overscreen the impu- 
rity spin. Additional conduction electrons must come in to remedy the over-performance. 
The process keeps going on resulting in a critical system. It is important to realize that the 
system is only critical along the time axis. As short time details are averaged out, the system 
approaches a universal limiting behavior given by the infrared fixed point. The impurity 
spin is asymptotically screened in the ground state which is a spin singlet. However, the 
infrared fixed point can be reached only in an infinite system in which the energy levels form 
a continuum. In finite systems, the typical distance between discrete energy levels, of order 
1/L for a system of linear dimension L, cuts off the scaling toward the fixed point. The 
ground state of a finite system has a residual spin and is doubly degenerated. A common 
definition of the residual spin is ^Tx(T) as the temperature T — > 0, where x(T) is the 
magnetic susceptibility. As we shall see, it is nonzero for a finite system. We shall also show 
that the entropy as given by the coefficient of the linear T term in the free energy is indeed 
\n(2S + 1) for a finite system. 

An intriguing feature of the multichannel Kondo model is that the entropy of an infinite 
system reduces to a smaller universal value [?,?,?] although the ground state is a spin 
singlet. If the entropy is still defined as the logarithm of the ground state degeneracy, the 
ground state could have a universal non-integer degeneracy. Since the entropy is deduced 
by calculating the free energy at small but finite temperature T which cuts off the scaling 
toward the fixed point, it is natural to relate the entropy to the effective residual spin at T. 
A necessary implication of this relation is that a tiny magnetic field will lift the degeneracy 
at T — > as we shall see later. 

Another interesting issue pertaining to all kind of Kondo problems, either overscreened or 
exactly screened, is the spatial size of the conduction electron screening cloud. This is a point 
underlying the resonant level nature of the Kondo problems. Considering the one-channel 
Kondo problem, it is well known that below the Kondo temperature T K the effective Kondo 
interaction enters the strong coupling regime. The energy gain from screening the impurity 
spin is of order Tk- If the screening cloud formed a localized bound state with the impurity 
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spin, it would have a spatial spreading of vf/Tk 3> 1/kF [?]• NMR experiment ruled out 
any conduction electron screening cloud bigger than one lattice spacing [?]. Although the 
physical argument for the screening cloud to have a size of 1/&f has been given [?,?], here 
we calculate the Knight shift for the multichannel Kondo model and explicitly show that 
the only length scale is 1/kp. We also explain why the same conclusion can be extended to 
the exactly screened case(k=l). 

The paper is organized as follows. In section [Q], the multichannel Hamiltonian and the 
Popov method are briefly recaptured. In section |T|, the general integral expression for the 
conduction electron self-energy to the order 0(k~ A ) is derived which can serve as a future 
reference. In section |V[ the integrals in the self-energy are evaluated at T = and the 
results are used to derive the RG equation and running coupling constant. In section [V], 
the scaling solutions for the conduction electron scattering rate and resistivity are obtained. 



The free energy is calculated in section |V| from which the specific heat and entropy are 
deduced. The magnetic susceptibility and field dependent magnetization are calculated in 
section |V11] for an equal spin gyromagnetic ratio for the conduction electron and impurity 



spin. The dynamical susceptibility for the impurity spin is calculated in section |V111| and the 



relaxation rate is deduced. The general case with different gyromagnetic ratios is considered 



in section [TX|. The Knight shift in the space surrounding the impurity spin is calculated in 
section |X|. The last section is devoted to a discussion of related issues. Part of the results 
has been briefly reported in [?]. Some details for the dynamical spin correlation functions 
are included in the appendix B. 

II. POPOV TECHNIQUE 

Without losing generality, we consider the local impurity to be a spin S = 1/2. The 
multichannel Kondo Hamiltonian in the magnetic field is 

fc,M=± k,k' f'T^ 

X=l-k X=l-k 



where S is a spin-1/2 operator, are the Pauli matrices, J is the Kondo interaction 
strength and Af is the number of lattice sites. We have set the conduction electron gy- 
romagnetic ratio and Bohr magneton equal to one so that the magnetic field h has the 
dimension of energy. We have introduced a parameter /is to account for possible difference 
between the conduction electron and impurity spin: us is equal to the gyromagnetic ratio 
of the impurity spin divided by that of the conduction electron. We shall adopt the usual 
cutoff scheme for the conduction electron band: —D < < D, with a constant density of 
states p per spin per channel. 

Representing the impurity spin in term of pseudofermions 

& = 7} E l]i°^vfv, 

fl,V=± 

Popov made an observation that the Kondo Hamiltonian without imposing constraint on 
the pseudofermions is disconnected in the pseudofermion charge sectors. That is 

Z = Tre-P H = Z + Z 1 + Z 2 = Ti5(n f )e' pH + Ti5(n f - \)e~ m + Tr6{n f - 2)e~ f3H . 
Among these three separate contributions, only Z 1 from the subspace 

n f = E ftU = 1 

a=± 

is physical. Popov's technique [?] is to add an imaginary chemical potential, iu = in/ {2(3), 
to the pseudofermions so that Z + Z 2 — 0. Using this trick, the partition function in the 
path integral formalism can be represented as 



Z 



= Jv[c,cJJ]exp -J Q P dT(Co + H) 



(2) 



£ o = E %A% A + E U( d r + iuo)U- (3) 

There is no additional constraint. The standard perturbation method then follows from 
this path integral representation. The impurities are assumed to be randomly distributed 
in space. The averaging over the impurity distribution is done according to the standard 
recipe as in the case of the spinless impurities [?,?] and we only keep contributions linear 
in the impurity density rii. The Feynman rules for constructing diagrams are listed in the 
appendix A. 



III. CONDUCTION ELECTRON SELF-ENERGY 



The first thing we shall calculate is the conduction electron self-energy in the absence of 
the external magnetic field. Up to the order J 4 , J 5 k and J 6 k 2 , the relevant diagrams are 
given in Figure |l|. The two subscripts in each term of the self-energy expansion indicate the 
powers of J and k respectively. After lengthy algebra, the final results at finite temperature 
are 



Efa n ,T) = E( 2)0 ) + £(3,0) + S(4,i) + £(4,o) + £(5,1) + £(6,2), 
3n,- „ r de 



£(2,0) (^n) = ~^ j2 P J 



lUJr 



tanh 



3ti f 

£(3,0) (^ n ) = -r^ j3 P 2 dedet 

16 J (iuj n - e)(e - e X) 



£(4,1) fan) 



3n. 



J (4,0) 



J 4 p 3 k J 

-jVj 



de de l de 2 f e €1( 



e - ex) fa n - e 2 )fa n - e 2 + e - e x ) ' 



<ie c?ei de 2 



, W n) = ^7T J P / /• 

32 7 (Kj n -e)(e-ei)(e-e 2 ) 



(2 



(4) 
(5) 

(6) 
(7) 
(8) 



'(5,1) l^n 



3n 
32 r 

+ ^J 5 A 

Id 



dex de 2 de 3 dei 



£2^3 



ei - e 2 )(e 4 - e 3 )(ei - e 2 + e 4 - e 3 )(iw n - e 4 ) 
dei de 2 de 3 de 4 tanh f <£>g ei 
(ei - e 2 )(e 4 - e 3 )(e 1 - e 2 + e 4 - e 3 )(z^ n - e 4 + e 2 - e x ) 

1 1 



de\ de 2 <ie 3 cfej tanh (^f 1 ) 



21n, 



128 



ei - e 2 )(iu; n - e 4 )(iu; n - e 4 + e 2 - ei) Ve 4 - e 3 e 2 - e 3 

dei <^ e 2 ^3 de^ tanh (^p) ^Sei 
(ei - e 2 )(e a - e 3 )(e 2 - e 3 )(zu; n - e 4 + e 2 - ei) ' 



(9) 



'(6,2) (^nj 



dei de 2 c?e 3 c?e4 des 



x 



X 



- — J 6 p 5 A; 2 
32 J fa n - e 5 )(ei - e 2 )(e 3 - e 4 ) 

(ia> n - e 5 + e 3 - e 4 )(ei - e 2 + e 3 - e 4 ) iu; n - e 5 + e x - e 2 + e 3 - e 4 
2 1 



£1^3 



+ 



(10) 



,«u; n - e 5 + e 3 - e 4 61-62 + £3-^4- 
where the integration range is always [-D, D]. Whenever necessary, the above integrals(and 
integrals throughout this paper) take principal value. We have introduced following short- 
hand notations, 
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^ = /(ei)/(-6 2 ) + /(-e 1 )/(e 2 ), 

= /(cs) /(-C4) /(-<*) + /(-cs) /(e*) /(e 5 ), 

= /(*0 /("<*) /(-cs) - /(-ci) /(c 2 ) /(en), 
= /(ei)/(-6 2 )/(e 3 )/(-e 4 )/(-6 5 ) + /(- ei )/(e 2 )/(-e 3 )/(e 4 )/(e 5 ), 

where /(e) = l/(e^ e + 1), is the Fermi-Dirac function. 

IV. RENORMALIZATION GROUP EQUATION AND SOLUTION 

Under the analytic continuation iu n — > + i0 + , £(u; + z0 + ) = S'(a;) + i£"(u;). The 
imaginary part of the self-energy £"(u;) is proportional to the total scattering rate of the 
conduction electron of energy u. At T = 0, the integrations in the imaginary part of the 
self-energy (p|)-([T0|) can be carried out and we obtain 

37T77/ 

D, g) = -^V (p + P x g In u + P 2 g 2 k In u + P 3 g 2 k + P 4 g 2 In 2 u + P 5 g 3 k In 2 u 
lop v 

+P 6 g 3 k\nu + P 7 g 3 k + P % g A k 2 \n 2 u + P 9 g A k 2 \nu + P w g A k 2 ) , (11) 

where g = Jp and uj = uj/D. All coefficients Pq to Pio are known: P — 1, P\ — —2, P 2 — 1 5 
P 3 = ln2-l, P 4 = 3, P 5 = -7/2, P 6 = 5-31n2, P 8 = 1, P 9 = 2 In 2 -5/2. The coefficients 
P 7 and Pio are not needed for deriving the RG equation to the sub-leading order. But they 
can be found from (^) and flTUp. Note that the g A order contributions from (|S|) contain only 
a ln 2 (d;) term. Since we expect the fixed point g* to be of order 1/k [?], our result ( [TTj ) 
includes all contributions up to the order 0{k~ A ). 

Since the dimensionless scattering rate ([□]) must be invariant under RG transformation, 
we obtain the following equation, 

[aL + ^S)^ D ^ = °- (12 > 

The beta-function, (3(g), of the Kondo interaction g can be deduced from (|ITD using the stan- 
dard RG technique [?]. The equation ( ]T2"| ) can be regarded as an equation which generates 
logarithmic series with fixed g, 
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Z"(co,D,g) = -(3(g) — J d(\nD)E"(u, D,g) + const. (13) 

Substituting ( pTTj) into (|13D, the integration over lnP> can be carried out easily. Expressing 
/3(g) to the sub-leading order as 

P(g) = 9 2 (*>i + K 29k + k 3 Q + ^g 2 k + n 5 g 3 k 2 ) (14) 

and substituting it into (|i~3"D , we obtain following equations for the coefficients Ki to k 5 from 
(|l3l) by equating order by order the coefficients of the polynomials on the left and right 
sides, 



Ki = 


Pl =-1, 

2P 


(15) 


K 2 = 


P 2 1 
2P 2' 


(16) 


K 3 = 


0, 


(17) 


K4 = 


P 6 -4P 3 «i l + ln2 


(18) 


2P 2 ' 


K 5 = 


P 9 - AP 3 K 2 1 


(19) 


2P 4' 



There are three consistency equations 

P 4 = -P1K1, 

P 5 = 2P 2Kl + ^P lK2 , 
P 8 = 2P 2 «: 2 , 

which are all satisfied. 

The obtained beta-function can be written explicitly as 

P{g) = -g 2 + k - 9 "+ k -{l + \n2)g*-^g\ (20) 

The intermediate fixed point g*, determined by /3(g*) = 0, and the slope of the (3(g) at the 
fixed point are then easily found to be 

fAh-*-¥\ Pi) 



k \ k 

r 

~k v ' /,• 



We find indeed g* ~ 1/k thus allowing a reliable expansion in the k — > oo limit. The slope 
which determines the critical exponent is universal. The perturbative result ( p2p agrees with 
the conformal field result 2/(k + 2) [?] up to the sub- leading order in the 1/k expansion. 
The running coupling constant Qr{uj) is determined by the differential equation 

^■ = «. (23) 
a In a; 

with the initial condition gn(u = D) = g. With our perturbative (3(g) of (ffiH), the solution 
of (|23|) covering the full range of energy scale from u <C T# to a; 3> Tr- can be obtained. 
We rewrite (|23|) in the integral form, 



a In a; = / — — - ~ g 



Id h M * J a (g') 2 (g'-g 

This leads to the full solution 

A 



l + g'ln2-^(g'r 



\gnH - g*\ = \g- g*\ [^) [ 9R (co)} kA/2 e~^\ (24) 
where Tk = Dg k ^ 2 exp(—l/ g). At uj < T Kl it has an asymptotic form of, 

9r{") = g*-c(^) A , C = (g*-g) (g*) kA/2 e~ A/9 * . (25) 

For an initial condition of weak coupling, g — > and D — > oo, ( = (g*) 1+/cA / 2 e~ A ^ g * . From 
(E3j), it is interesting to note that the running coupling constant has a power law behavior 
not only at low energy uj < Tk but also at high energy uj > Tk, underlying the critical 
nature of the system. 



V. SCATTERING RATE AND RESISTIVITY 

To obtain the scaling solution for the conduction electron scattering rate pS", we can 
use the RG invariance, 

D, g R (D) = g) = D', g R (D')). (26) 

Choosing D' = u>, the logarithmic terms in pTl"(u>, u, gnioj)) of flO] ) drop out and the scaling 
form for (|TTD is 
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pE"(uj,D,g) 



37rrii 



gl(u;)-(l-\n2)kg 4 R (u) 



37m,- 



4(k + 2)< 



1-Ar 



CO 



T, 



K 



(27) 



where A = k([l — (4 — 6 In 2) /A;]. A similar RG procedure will be performed repeatedly 
later. Note that instead of a Lorentzian frequency dependence in the exactly screened Kondo 
problem, pS" ~ T^/iuj 2 + T|-), the scattering rate ([27]) has a cusp at ai = 0. At nonzero 
temperature, the cusp is expected to be smoothed out. 

In order to calculate the resistivity, we need the temperature dependent imaginary part of 
the conduction electron self-energy S"(u;,T). Identifying the transport relaxation time due 
to the Kondo exchange as r ex (u, T) = l/[2H"(u, T)] since there is only s-wave scattering [?], 
the total scattering rate is 1/t(u, T) = 1/tq + 1/t cx (oj, T), where To is the ordinary relaxation 
time in the absence of the Kondo interaction. The ordinary scattering 1/tq could arise from 
spinless impurities or defects which are assumed to be located at different lattice sites and 
uncorrelated to the impurity spins. A low impurity spin density is assumed, rij <C 1, such 
that r <C r ex . The total relaxation time is substituted into the conductivity expression [?], 

t(w) n e e 2 f°° , t ( A t 



a(T) 



n.e* 



du- 



ctal- 



(28) 



m e Jo 2Tcosh 2 (^) m e Jo 2Tcosh 2 (^) \ r cx (aj) J ' 
with e, m e and n e denoting the conduction electron charge, mass and density respectively. 
With our perturbative expression for H"(oj,T) to the order k~ 3 , we find for the resistivity 
due to Kondo scattering, 



37T m P 



00 dx 



16 n e e 2 p Jo cosh % 
3tt m e 



d/t tanhf - T / D 

l-g dy l-g 2 k^(x,- 

d/t x — y V 1 



8 n e e 2 p 



rug 



1 + 0(g,g 2 k) + 0(g,g 2 k)x In 



D 



(29) 



where 



D/T I'D IT — l - i 

dy 1 dy x ^^ 
-d/t J -d/t [y - yx) 1 



1 



1 



1 



(30) 



e yi-y-x _|_ I e yi-y+x _|_ I 

Some difficult integrals have to be evaluated to obtain the sub-leading terms in fl29D . So we 
limit ourselves to the leading order. Performing the same RG procedure as (|27|), we obtain 

37T 2 71,- 



4k 2 n, 



■P. 



(o) 



rp x A 



(31) 
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where pf^ = Ait / (e 2 kp), is the resistivity in the unitary limit and kp denotes the Fermi wave 
vector which is related to the conduction electron density of states through p = kpm e / {2it 2 ). 
The T = value of the resistivity and the exponent A have been reported [?] and an exact 
expression for the resistivity for k = 2 up to a constant factor ( has been derived recently 
[?]. As in the exactly screened Kondo problem, the corresponding resistivity decreases upon 
increasing temperature. 

An important point is that even at T = 0, there is still inelastic scattering. In other 
words, the impurity spin together with its asymptotical screening cloud(see magnetization 
section) cannot be viewed as inert. Let's see how the assumption of elastic scattering only 
at T = would run into trouble. Following [?], the total resistivity can be divided into 
elastic and inelastic contributions, 

dp e = dp? + cos(20 o ) dp 1 :, (32) 

where <fro is the scattering phase shift at the Fermi energy. This separation is valid at least for 
weak inelastic scattering which would indeed be the case at low T if only elastic scattering 
were present at T = 0. For elastic scattering, we could identify 

pS"(cu = 0,T = 0) ~ n i S ( S + 1 * ) sin 2 0o, 

It 

which would lead to 

% ~ ~- ( 33 ) 



k 



The resistivity due to the elastic scattering can be calculated by using (p7|) as the scattering 
rate. To the leading order, 

du 3it 2 rij 







(34) 



Substituting (|3l|), (p3|) and (|3"4]) into (|32"D, we would obtain the inelastic contribution to the 



resistivity dp 1 : < which certainly is impossible. This contradiction indicates that there are 
both elastic and inelastic scattering at T = 0. For the two channel case, k = 2, it has been 
suggested that there is actually only inelastic scattering [?]. 
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VI. SPECIFIC HEAT AND ENTROPY 



To calculate the free energy, we can use the linked cluster theorem, 

oo 

F = F - U n , 



(35) 



n=2 



where F is the free energy of the non-interacting Fermi sea. The diagrams for the first 
three U n are drawn in Figure [2[ The results after completing the Matsubara frequency 
summations are 



16 y J e-ei 



U 3 = k I de ~t de 2 de 3 
3n, ; 



/(ei)/(e 2 )/(-e 3 ) 
(e 3 - ei)(e 3 - e 2 ) 



+ 2 Permutations 



Ui = ^9 4 k 2 J dei de 2 de 3 de A /(-ei) /(e 2 ) /(e 3 ) /(-e 4 ) 

'j _ g/3(e2-ei+e3-e4) g/3(c2-ei) _ g/?(e3— 64) ' 
+ 



(36) 
(37) 

(38) 



(ex - e 2 )(e 3 - e 4 ) \ e 1 - e 2 - e 3 + e 4 e x - e 2 + e 3 - e 4 < 

As noted by Kondo [?], it is quite delicate to extract the linear T terms in ( ]3"T| ) and 
Here we follow the method used by Kondo and present the calculations only for the integrals 
not calculated in [?]. 

The calculation of £/ 2 is simple and straightforward. The result is a constant plus T 2 
corrections. To calculate U 3 , we define a new integral J 3 , 

f(ei)f(e 2 )f(-e 3 ) 



D 



dex de 2 de z - r— — , 

d (e 3 - ex)s (e 3 - e 2 )s 



where we have introduced the notation 



1 



x 



(x)s x 2 + 5 2 

The 5^0 limit is taken after the integrations are completed. Kondo showed [?] 

3m 



U* 



7T 



(39) 



so we only sketch the calculation for I 3 . 
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*3 [ . ^ ^ /(6x) /(6 2 ) /(6 3 ) 



f deide 2 f(e 1 ) f(e 2 ) [ -, r-^ r - / de 1 de 2 de 3 - 

J J-d es — ei )x e* — eg ■ * j i 



'-£) (e 3 - ei){ (e 3 - e 2 )<5 i (e 3 - ei) s (e 3 - e 2 )a 

The first integral has no linear T contribution. Thus, up to a constant term, we have 
1 



1 r u 

h = —z\ de x de 2 de 3 /(ei) /(e 2 ) /(e 3 ) 

3 J-D 



1 1 1 

+ 7 ; r- + 



(e 3 - £1)5 (e3 - £2)5 (ei - e 2 ) 5 (ei - e 3 ) 5 (e 2 - ei) 5 (e 2 - e 3 ) 5 
. , f(e 1 )f(e 2 )f(e 3 )[(e 1 -e 2 ) 2 + (e 1 -e 3 )(e 2 -e 3 )}5' 



/ dei de 2 a!e 3 

J-D 



3 J— Z) [( ei _ e2 )2 + 5 2 ][(ei _ e3 )2 + 5 2 ][(e2 _ e3)2 + 52] 

7T 2 

= — 5- / dei de 2 de 3 /(ei) /(e 2 ) /(e 3 ) 5(ei - e 3 ) 5(e 2 - e 3 ) 

3 J-D 

71-2 .£> 2 

Combining the above result with (|3~9| ), we obtain 



r^vr 2 ., 



C/3 = const. ^-0 * T - (40) 

Next, let's turn to U4 which can be calculated in a similar way. Defining a new function 



h = J dei de 2 de 3 de 4 /(-e x ) /(e 2 ) /(e 3 ) /(-e 4 ) 



J / 1 — e/ 3 ( e2_<E1 + <E 3 _e4 ) g/3(«— ei) _ g/3(e3-e4) 

' + 



(ei - e 2 )5(e3 - £4)5 \ (ei ~ e 2 - 63 + £4)5 (ei - e 2 + £3 - £4)5 

f(-ex)f(e 2 )f(e 3 )f(-e A ) 



4 y dei ^e 2 de 3 de± 



(ei - e 2 ) 5 (e 3 - 64)5(61 - e 2 - e 3 + 64)5' 
then following the same steps Kondo employed to prove (|5§D, it's straightforward to show 

Ua _ ^± g ±k 2 h = ^-^g 4 k 2 T. (41) 
64 y 32 y v ; 

It can be further verified that J4 has no contribution linear in T. Therefore, 

U A = const. + ^Vfc 2 T. (42) 
32 



Substituting (fiO|) and ( |4"2] ) into (p5[) and setting rij = 1 for simplicity, we obtain the free 
energy shift due to the presence of the impurity spin, 

F imp (T) = -E - Tin 2 + (V - jj*y) + C% 4 fc), (43) 

13 



where E is the ground state energy. 

Since the free energy shift is RG invariant up to an additive constant, we obtain the 
scaling solution for the free energy shift at T — > after carrying out standard RG procedure, 



■2 



F imp (T) = -E' -Tln2 + ^T 



k9%(T) - \k 2 g%{T) 



= -E' -T{\n2-—)- — C 2 T + 0[T(T/T K ) 3A ] , (44) 

where E' is in general different from Eq of (f|3). The impurity specific heat is 

r T 2 F imp 3vr 2 2 /Tf 

Cimp = - r ^ = T (A fej ■ (45) 

The critical exponent, a = 2 A, agrees with the previous result [?,?]. The impurity entropy 
is reduced to a universal value [?] 

5W(r = 0)=ln2-^~lB^2-^. (46) 

Actually, the correction to the entropy occurs only when one first takes the thermodynamic 
limit M —>■ oo and then the limit T —>■ 0. For a finite system, the integrals in {7 3 and 
should be replaced by discrete momentum summations over the Brillouin zone which will 
not give rise to contributions linear in T to the free energy, thus no correction to the bare 
In 2 term of the entropy. Therefore, the finite system always remains doubly degenerate. 



VII. FIELD DEPENDENT MAGNETIZATION 



A natural question following the finite entropy is whether or not there is a corresponding 
residual impurity spin. To answer this question we calculate the field dependent impurity 
magnetization defined as the total magnetization of the system subtracting the free Fermi 
sea contribution. In the presence of the external magnetic field, the leading order diagrams 
to the free energy are shown in Figure |3|. This is an infinite series. The diagram with I 
vertices has a contribution, 



SFi 



1 k l 

'p7 { 



-J 



Af(3 



EE 



-a/2 



CT=± iuj n -e % -ah 



E 



-a'/2 



uJ m ,o-'=± 



iuj m — icuo — a'ush 



(47) 
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The momentum summation over the conduction band is readily carried out, 
J 



a 



A/7? £ iu n - e % - ah 



g [ D de[f(e + h)-f(e-h)] = -2gh. 

J-D 



(48) 



Summing up all the diagrams in Figure |3] and using (fi8|) , the free energy shift due to the 
impurity spin is 

oo 

F imp (T,h) = '£6F l = 

l <T— ± ( = 1 \ " / L J Z = iLOo+<Tflgh 

Introducing the primitive of f(z): du(z)/dz = f(z), the series can be summed up, 
F imp (T, h) = [u{iu + (Tfi s h - akgh/2) - u(iu + a^ s h)] . 





W-'fiz)] 


^—agkh^j 1 


dz 1 - 1 



The corresponding magnetization is 



5M(T,h) = -^ 



(3h 



( ^ 



— tanh(Pnsh). 



Adding to the above result the bare magnetization of the impurity spin which just cancels 
the second term of the above expression, we finally obtain 



M(T,h) = Ls- yjtanh 



Ph\ns- 



kg" 



(49) 



This is the leading order temperature and field dependent magnetization. 

In this section, we shall mainly consider the case us = 1, i-e. an equal gyromagnetic 
ratio. The expression ( |49| ) has the form characteristic of a reduced free spin, 

M = S cS ttmh(S c sh/T). 

The fact that S e ff —>■ following a power law implies that the impurity spin is completely 
screened only for the infinite system. A finite system of size L behaves as if there is a 
partially screened spin S e fi ~ L~ A since there is a minimum energy unit 1/L to cut off 
the scaling. It also suggests that the ground state degeneracy will be lifted by any small 
magnetic field at T = 0. Using the Maxwell relation dS/dh = dM/dT, the leading order 
magnetic field dependence of the entropy is 
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S imp (T,h) =ln2+ / 
Jo 



•> dh ,8M<T,H) 



dT 



f f3h\l~gk/2\ x IT -> 

= ln2- / — -« >0, for ^ (50) 

Jo cosh x ^ / g 

We expect that this is true to any order. In particular, the entropy change, 7r 2 /(2k 2 ), in the 
next order (|46|) will be removed by the magnetic field. 

From fl49|), we obtain the low temperature magnetic susceptibility and its scaling form 



kgV a (, kg R (T)\ 2 T<TK fkCY 1 



indicating that the fixed point is a spin singlet since S 2 S ~ Tx imp (T) ~ T 2A — > at 
T —>■ 0. The spin is quenched very slowly compared with the linear dependence in the exactly 
screened case. Also from flj9p, we determine the field dependence of the zero temperature 
magnetization [?], 

MiT = 0, h) = 1 - f = 1 - h -if> f (A)" . (52) 

This shows that the spin operator has the scaling dimension A. 

Using the well known results for the bulk specific heat and magnetic susceptibility: 
Cbuik = 2krr 2 Tp/3 and Xtmik = 2fcp respectively [?], the leading order Wilson ratio is deter- 
mined from and (51), 



W = Ximp Cbulk = (f) ? fe) |fe7T 2 Tp = fc 3 

CimpXbulk 3|£ C 2 A ^) 2A 2fc P 36' ^ ' 

in agreement with the conformal field result [?]. The Wilson ratio is universal because the 
only parameter having a possible dependence on the cut-off scheme is ( which is canceled 
out. 

For lis ^ 1, the external magnetic field couples to an unconserved spin operator. The 
magnetization will acquire an anomalous dimension and we shall discuss it in detail in 



section IX 
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VIII. DYNAMICAL SUSCEPTIBILITY 



For the one-channel Kondo problem, it is well known that the impurity spin flipps at a 
typical rate of Kondo temperature Tk- At T < Tk, the impurity spin is effectively inert 
leading to a fixed point of local Fermi liquid type described by a phase shift. To understand 
the multichannel case, we calculate the dynamical spin susceptibility. As we shall see, the 
typical spin flipping rate is given by the temperature. 

We proceed as usual by first calculating Matsubara spin-spin correlation function. 

Xf (r) = (TS(r) • 5(0)> = ± £ Xf(™n) , (54) 

where v n = 2rwr//3. To the sub-leading order, the diagrams are shown in Figure f|. The 
result of calculating these diagrams is 

X/K) = <W>-r I 1 ~ 9 2 H^PD - In 2 + h)\ - -M—K{iv n \ (55) 



4 

where I = / °° rfx[tanh(x/2)/x - 1/(1 + x)} ~ 0.125 and 

rD rD de 1 de 2 



K{iv n 



-dJ-d [e t - e 2 y + K 



f^) - f(e 2 ) , 

- J {£l)Onfl 



si - e 2 
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;i-5 n , ) + O( J D- 1 ). (56) 



To understand the result we just derived, let's consider x{ z ) i n the complex frequency plane 
iv n — > z outside a disk of radius of order T (Figure [5]). Knowing its values at a set of discrete 
points on the imaginary axis, z = iu n , is enough to give a unique analytic continuation, 

3t 2 1 3?r 2 % sgn(Imz) 

X/Fn = -r£ fe, — r -> -rff *: , 57 

4 |v n | 4 z 

where sgn(Im z) means taking the sigh of the imaginary part of z. When z approaches the 

real axis, z = u + id, Xf{ u + i$) = X/( w ) + * sg n (<5) x"f{ u )i an d we find 

XJ( W ) = ^s 2 Jfc 1. (58) 

This is valid in the whole plane outside a disk of radius of order T. There is no energy 
scale to cut off 1/u dependence. Certainly for w<T, this dependence should be flattened 
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out as seen from the v n = expression of fl55|). If we recall that in the exactly screened 
Kondo problem(/c = 1), }("j{uj)/uj ~ l/(u 2 + T K ), the system gives sizable magnetic response 
only when all the contributions up to the energy scale Tk are included. In other words, the 
impurity spin is effectively quenched on time scale longer than 1/Tk- In the multichannel 
case, the role of Tk is played by the temperature T! The impurity is only marginally 
quenched since the time scale is 1/T. 

In the other limit uo < T, carrying out the analytic continuation as in [?], we find 

X '^ T) ~ T 2A_1 . (59) 

Note that if us ^> 1, the coupling of the nuclear magnetic moment to the impurity spin 
dominates and the above result is proportional to the NMR relaxation rate 1/(T\T). In 
the limit k ^> 1, local probe basically sees a nearly free impurity spin. In the two-channel 
case, k = 2, 2A = 1, we see that flSSP will have no power law divergence in T except possible 
logarithmic dependence. 

In the two-channel case, numerical work in the non-crossing approximation [?] and the 
solution at a special Toulouse point [?] have found x'f( u ) ~ tanh(u;/2T)/(u; 2 + T K ). The 
tanh(co>/2T) term has also been produced by the conformal field theory method and is a 
property of the fixed point itself. It would be interesting to see if the remaining Lorentzian 
form with width T% can be reproduced from the perturbation of the leading irrelevant op- 
erator in the exact conformal field theory calculations. For k > 2, we do not expect a 
Lorentzian form with finite width Tk for the following reason. The multiplicative renor- 
malization factor for X/( see next section) will bring an anomalous dimension u 2A in (|58|). 
We see that ~ w 2A_1 for T < u < T K . Since 2A = A/{k + 2) < 1, this frequency 

dependence cannot be reconciled with a Lorentzian form with a finite width. 

IX. MAGNETIC SUSCEPTIBILITY FOR fx s ^ 1 

Besides the time ordered impurity spin-spin correlation function studied in the last sec- 
tion, we can also define 
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X fe(q,r) = (TS e (q,r) • S(r = 0)) = -£ X/e (g> n ) ^ 

P n 

Xe (q,r) = (TS e (q,r) • £(0,0)) = ~5> e (S>n) e"^ T , 

P n 

where the conduction electron spin operator is defined as 



(60) 
(61) 



k \=l-k 



(62) 



The Feynman diagrams for these two correlation functions are shown in Figure |] and cal- 
culated in the appendix B. 

When the impurity spin has a different gyromagnetic ratio from the conduction electrons, 
a uniform external magnetic field couples to an unconserved spin operator Sh = S e (q = 
0) + UsS. The magnetic susceptibility is then only RG invariant up to a multiplicative 
renormalization factor, 



Ximp (T,g,D) = [Z h (g,g R (D'))} 2 x imp (T,g R (D'),D 



(63) 



where Zh is the multiplicative renormalization factor for the operator Sh and we recall 
gu(D) —g. If fi s = 1, then Z h = 1. The RG equation for \ is 



Xim P {T,g,D) = 0, 



ih(g) = d 



dhiZh 
dD 



(64) 
(65) 



The perturbative result for \ to the sub-leading order can be found from the results for 
three dynamical correlation functions calculated in the last section and appendix B, 

4 r . 



Xim P (T,g,D) 



H Us 



Xe[q = 0, U n = 0) + 2/J, S Xfe(q = 0,Vn = 0) +/i S X/(^n = 0) 
r \ 2 / 



1 - g 2 k(\n /3L> - In 2 + J ) - g 2 k fi 8 



gk 



In 2 } . 



(66) 



The leading term of jh{g) can be obtained from ( |66"D using the same method of section [TV] , 



ih{g) 



fIS ~ 1 g 2 k. 



tys - 9 k 



(67) 
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Its value at the fixed point, jhig*), gives the anomalous dimension of the operator Sh- The 
anomalous dimensions for the operators S and S e (q = 0) are the limiting values at fj,$ — > oo 
and us — > respectively. 

To obtain the scaling solution, we first find the multiplicative renormalization factor, 

7i m<n r m ,M9') ( gk-tys Yf 9R{D')k-2 V 

Setting D' = T in ( |63|) and (|68[) , we obtain the leading order susceptibility(since we only 
have leading order Z h ), 

Ximp (T, g, D) = {^zt) 2 [Xtop(T, ^ ^)], s =i • (69) 

This is the main result of this section: the static magnetic response is the same, up to a 
constant factor, no matter whether or not the conduction electrons and the impurity spin 
have the same gyromagnetic ratio. I believe this conclusion remains true for all overscreened 
cases and even for the exactly screened case of k — 1. If the initial starting point belongs to 
weak coupling, g — > 0, 

Xim P (T, g, D) = /4 [Xim P (T, g, £>)] Ms=1 • (70) 

This is the well known result for all kind of Kondo problems: the contribution to the 
magnetization from the conduction electrons is suppressed by a factor T^/D compared with 
that of the impurity spin [?]. Extending the treatment in this section to the field dependent 
magnetization and dynamical susceptibilities is straightforward. 

X. KNIGHT SHIFT 

Knight shift gives direct measurement of the spatial structure of the conduction electron 
screening cloud. Consider an impurity spin sitting at the origin and a uniform magnetic 
field h being applied to the system, Knight shift measures the magnetization in the space 
surrounding the impurity spin, 
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M{f)= £ a « A (rO c CT)A (r)) = X (r) h, xW^E^I^ (71) 



<7 = ± 

A=l-fc 



We implicitly assume that the free Fermi sea contribution to M(r) has been subtracted. The 
result for the Knight shift is contained in the two general dynamical susceptibilities defined 
in the last section, 



Xio) = 3 [Xe(5,i^n = 0) + iisXfe(q,iv n = 0)] ■ 
Using the results from the appendix B, we find to the sub-leading order in 1/k, 



(72) 



x(q)--^MQ)\Us- 9 Y 



l-g 2 k(ln^-\n2 + I ) 



Mo) = 

ni(g) = 



6 k 6 k+q 



tanh || [/(eg) - /(e^) /(e^>) - /(e^ + -) 



1 1 



(73) 
(74) 

(75) 



where J ~ 0.125, given in the last section. We have included proper normalization factors 
in the definitions of n (<7) and Hi(q) such that U o (q = 0) = 1 and U^q — O) = 1. Note that 
the Fourier transform of Ho(q) has the well known Friedel oscillation form in space [?]. The 
Fourier transform of Hi(q) should give similar spatial variation. The key point is that their 
Fourier transforms only depend on kpr, where r is the distance to the impurity spin. 
Since x(o) satisfies the RG equation, 



D dD + ^ 9) dg + lh{9) + frfc($)U=o 



x(q) = o. 

The scaling solution can be found in a similar way to the last section, 



(76) 



X(q, T, g, D) = Z h (g, g R {T)) [Z h (g, g R (T))]^ =0 X (q, T, g R (T),T) 



(77) 



Using the results for Zh and g R , we obtain to the leading order, 

X(q, T, 9, D) = Qk ^Zly S) X^{T) U (q), 

21 



(78) 



where Xim P (T) is the uniform susceptibility found in section [VI]] . A knowledge of Zh to the 
sub-leading order could extend ([T8|) to the next order. When Fourier transforming ( |78| ) to 
the real space, M(r) is a function of kpr. So the only length scale is Thinking over 

why the length scale vf/Tk does not appear in x(<l), we observe that n and 111 are functions 
of q/kp. A new length scale Vp/T^ would appear only if the logarithmic singularity were 
cut off by VFq in the limit vfQ > T. This doesn't happen! It means that the multichannel 
system is critical only in the time direction, not in the space. 

We can extend the statement that there is only one length scale 1/kp to the one-channel 
problem. In this case, the effective Kondo interaction flows to strong coupling at low energy. 
From the renormalization group point of view, it means that summing up leading(or leading 
plus sub-leading) logarithmic series is not enough. But reorganizing the perturbation expan- 
sion in the coupling constant g into successive logarithmic series(the first term of each series 
has successive higher power in g) according to the renormalization group is still possible. 
This is equivalent to write (3(g) as an expansion in g. The occurrence of a strong coupling 
fixed point means that we need all terms in the expansion of (3(g) to describe the low energy 
behavior. Suppose we could find all of them(infinite terms), we would have a correct low 
energy theory by summing up all terms. The fact that the infrared singularity is not cut 
off by VFq does not change merely because we have to include high order logarithmic series. 
The necessary consequence is that Vp/Tx will never have the chance to appear as a length 
scale. Certainly, if the flowing of an effective interaction to strong coupling leads to a phase 
transition like in higher dimensions, the above argument may not hold. But the impurity 
problem has a dimension + 1, prohibiting a phase transition. Thus, our calculation pro- 
vides an explicit analytical demonstration that the screening cloud has a spatial size of order 
1/kp, and agrees with the well known experimental result [?]. 
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XI. DISCUSSION 



We have carried out a comprehensive study on the non-Fermi liquid fixed point of the 
multichannel Kondo model emphasizing several intriguing aspects. As the actual realization 
of this model in heavy fermion and metallic glassy systems is still an open question, a 
thorough understanding of the model should be very helpful in devising new experimental 
tests and making comparison with experimental results. Part of the driving force behind 
the recent resurgence of interest is due to the resemblance of its low energy behavior to 
the normal state properties of high Tc cuprates. This similarity has been further exploited 
recently [?] arguing that a situation similar to the two-channel model is realized in the 
copper-oxide plane. 

One interesting followup subject is to study the possibility of converting the local non- 
Fermi liquid fixed point into a coherent lattice one by forming an overscreened Kondo lattice, 
especially in two or three dimensions. The first step to study the two-channel-two-impurity 
problem has been carried out using the numerical renormalization group method [?]. The 
one-impurity fixed point is found unstable against developing correlations between the im- 
purity spins. This is expected in the view of the asymptotical screening of the impurity 
spin. One artificial way to suppress magnetic correlations between impurity spins is to go 
to infinite dimension. Unfortunately, the presence of finite degeneracy on each site prevents 
development of true coherence. In any realistic situation, nontrivial magnetic correlations 
must intervene to lift the residual degeneracy of the impurity fixed point. If one wishes to 
follow the successful route of the heavy fermion theory from impurity to lattice once again, 
in the current intensive search of non-Fermi liquid models to describe the normal state of 
high Tc cuprates, only a non-degenerate fixed point of an impurity model has the chance to 
succeed. Development along this direction in search of a promising starting point of impurity 
model has been reported recently [?]. It is worthwhile mentioning that a quantum critical 
system with a Fermi surface has a true coherent non-Fermi liquid infrared fixed point in 
high dimensions. An example is an electron gas coupled to transverse gauge field where a 



23 



self-consistent solution of the infrared fixed point has been derived [?] . It is very interesting 
to note that this system is quantum critical at T = in 2D and 3D but without developing 
long range order. 
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APPENDIX A: FEYNMAN RULES 

To be self-contained, we list the rules for constructing Feynman diagrams in the presence 
of an external magnetic field. These rules also define our convention. 

1. For the contributions of the nth order in the Kondo interaction, 0(J n ), draw all 
topologically distinct diagrams with n vertices. 



2. For each conduction electron propagator, draw a solid line, 



k, iuj r , 



G$(k,icu n ) 



-5. 



(Al) 



— * 

There is a summation over each momentum k of the internal conduction electron 
propagator. For each pseudofermion propagator, draw a dashed line, 



k,iu n 



-5 



[IV 



iuj n — iu — nush 



3. Each vertex is associated with a factor, 



(A2) 




J 7} 



k, iu>„, u 



uv m ,v 
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(A3) 



4. Each independent internal frequency is summed over. 

5. Each conduction electron loop contributes a factor — k and each pseudofermion loop 
contributes a factor —1. 

6. For a diagram of order J n , there is a numerical factor: (Number of different 
connections) jn\. This combinatorial factor will be given explicitly in the figures of 
this paper for the diagrams we calculate. 



APPENDIX B: DYNAMICAL SPIN-SPIN CORRELATION FUNCTIONS 



Three dynamical correlation functions have been defined in the sections [VIII| and |X[ To 
the order 0(1) + 0(1/ k), their diagrams are shown in Figure |j. Each diagram contains two 
external vertices which are represented by two open ends at the left and right sides. The 
external vertex at each end is joined by two propagators. If the spin indices of these two 
propagators are <j\ and 02, the corresponding vertex is associated with a vector of matrices 
s a2U1 . The two vectors of matrices from the two external vertices at the two ends form a 
scalar product s a2(T1 ■ s Uia . A . In Xe, only one of the two external vertices has the two joining 
conduction electron propagators differing in their momenta by q. 

The result for Xf(i p n) is given in fl55|). The results for the other two are 

Xfe&ivn) = -<W ^gk/3{u (q) [l-g 2 k(\n(3D-\n2 + I )] + g\n2U 1 (q)} 

+(l-6 n fl)xfe(.<t,iv n ), (Bl) 

g k 3 In 2 

X e (q, iv n ) = — — 5 n ,o Xfe(q, iv n ) + —^-g 3 k 2 (3Tl (q) - (1 - 5 n ,o)Xfe(q, iv n ), (B2) 

where n and TEi are defined in section [X| and for u n 7^ 0, 

3 

Xfe(q,iv n ) = -g 2 k [gkUo(q)K(iv n ) + 2K(q,iv n ) - 2L(q,iv n )] , (B3) 

o 

^ k,k' k k'+q^k t k'J ^ u n 

L ^ ™n) = -^E , — )2 2 - (B5) 

^ k %' k k ' ^ k ' k'+q) ^ U n 
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Note that K(q — 0, iv n ) = K{iv^), is a generalization of K{iv^) to finite q. In the expressions 
for Xfe and Xe, we have dropped contributions of order 1/D to the u n = Q components with 
respect to j3. 
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FIGURES 

FIG. 1. Feynman diagrams for the conduction electron self-energy. 

FIG. 2. Diagrams for the free energy. The solid line represents the bare conduction electron 
propagator. The big circles stand for self-energies. 

FIG. 3. Diagrams for the free energy in the magnetic field. 

FIG. 4. Diagrams for three dynamical spin correlation functions. 

FIG. 5. Analytic continuation in the complex frequency plane carried out in the region outside 
a disk of radius of order T. 
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